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slicing condition will be useful for simulating higher-dimensional black-hole spacetimes in numerical 
relativity. For the case of D = 5, we present analytic solutions of the intrinsic metric, the extrinsic 
curvature, the lapse function, and the shift vector for the foliation by the stationary maximal 
hypersurfaces. This data will be useful for checking five-dimensional numerical relativity codes 
based on the moving puncture approach. 
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I. INTRODUCTION 



The merger of binary black holes is one of the most important sources of gravitational 
waves for gravitational-wave detectors. As a result of longterm effort by many numerical- 
relativity researchers, it is now feasible to theoretically predict the orbital evolution and 
merger process of astrophysical binary black holes and resulting gravitational waves emitted 
from these systems by numerical relativity. After the first report of longterm simulations of 
binary black holes by Pretorius several groups have also succeeded in longterm simula- 
tions . There are basically three formulations in numerical relativity, workable 
for simulating dynamical black-hole system; one is the "generalized harmonic" formulation 
with the help of "black hole excision" l|, second one is the so-called Baumgarte-Shapiro- 
Shibata-Nakaumra (BSSN) formulation [g, |9|] together with the "moving puncture" method, 
and third one is the hyperbolic formulation with the black hole excision 6|. 

Another scenario of black hole formation in the system of two relativistic objects was 
pointed out in the non-astrophysical context 1^, 
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12| : mini-black-hole production at 



high-energy particle collisions in particle accelerators in the framework of the brane world 
scenarios (i.e., the so-called TeV gravity scenarios) Q- Motivated by this possibility, 
mini black holes in the particle colliders have been studied by many researchers in the 
past decade from a wide viewpoints (see [isl for a recent review). There are two marked 
differences between the black hole production in the colliders and the black hole mergers 
in astrophysical situations. One is the relative speed of the objects. In the high-energy 
particle collisions at the LHC, the 7-factor of an incoming proton is ~ 7 x 10^, and that of 
the partons consisting of each proton can be much larger. The other is the number D of 
the spacetime dimensions, because the spacetime dimensions higher than four is essential 
for the TeV gravity scenarios. 

In order to understand well the black hole production in the particle colliders, numerical 
relativity is probably the unique approach. To perform a simulation of this system, two 
techniques have to be developed: One is the technique to handle the high-energy objects and 



the other to simulate higher-dimensional 
developed for the four-dimensional case 
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ack-ho 



17 



e spacetimes. The first technique has been 



18| by modeling the high-energy two-particle 



system as the high-velocity two-black-hole collision and extrapolating the results to the 
ultra-relativistic regime. Such simulations have not been performed for higher-dimensional 
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5. Q,""^ 



spacetimes yet (but see Refs. [19|, l20( for studies on the apparent horizon in the system of 
two high-energy particles in higher dimensions). Also a numerical code to simulate higher- 
dimensional spacetimes with the BSSN formalism has been recently developed 2l| (see a 
list of proposed several test simulations and the successful results). 

In this paper, we investigate the maximal slicing condition for the D-dimensional spher- 
ically symmetric black- hole spacetime {D > 5), known as the Schwarzschild-Tangherlini 



solution 



22j. This study is motivated by numerical relativity performed with the BSSN 



formalism and the puncture method, by which whole spacetime region including the black 
hole interior is simulated, and thus, the slicing condition that has the singularity avoidance 
nature has to be adopted. For the four- dimensional Schwarzschild black hole, it was shown 
that the sequence of maximal slices never plunges into the curvature singularity but asymp- 
totes to the so-called limit surface 23|]. We shall show that the sequence of the maximal 



slices avoids the singularity also in the higher- dimensional Schwarzschild-Tangherlini space- 
time. This gives the theoretical foundation that both the maximal slicing condition and the 
puncture gauge condition (see 24, l25|]) have the desired features in higher-dimensional nu- 
merical relativity. Another important aspect of this study is to derive a stationary sequence 
of the maximal slicing hypersurfaces that provides a useful analytic solution for a benchmark 



test of higher- dimensional numerical relativity codes, as discussed in Refs. 25|, |26| in the 
four-dimensional case. 

This paper is organized as follows. In Sec. II, we derive the general solution in the 
maximal slicing for the D-dimensional spherically symmetric black-hole spacetime. Then, 
in Sec. Ill, we investigate the foliation by the maximal hypersurfaces which are reflection 
symmetric at the wormhole slot. In Sec. IV, we study the foliation by the stationary maximal 
hypersurfaces. After explicitly describing the analytic solution in the case D = 5,we show its 
usefulness as a benchmark test for numerical relativity, by performing a numerical simulation 
adopting the analytic solution as the initial data. Sec. V is devoted to a summary. In 
Appendix A, the Kruskal extension of the Schwarzschild-Tangherlini spacetime is analyzed. 

We adopt the geometrized units c = G = 1 throughout this paper, where c is the speed of 
light and G is the gravitational constant of a D-dimensional spacetime. The Greek indices 
(/i, u, ...) represent the components of a spacetime, while the Latin indices ...) represent 
the components of a space. 
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II. GENERAL SPHERICALLY SYMMETRIC MAXIMAL SLICING 



Following Estabrook et al. [23|, we derive general spherically-symmetric maximal slicing 
of the Schwarzschild-Tangherlini spacetime, i.e., the D-dimensional spherically-symmetric 
vacuum black- hole solution with D > 5. The general form of its line element is 

ds^ = -a^df + -f{dr + -i-^pdtf + r^dQl_^, (1) 

where dQ'jj_2 is the line element of a (D — 2)-dimensional unit sphere with the {D — 2)-area 
■= {D — l)n^^~^^^'^ /T{^^j^), and a, /3, and 7 are functions of t and r. 

Einstein's equations with appropriate coordinate conditions lead to the basic equations for 
a, (3 and 7. Because the {D — 2) -dimensional spherical-polar coordinate system is uniquely 
determined in the spherically symmetric spacetime, two conditions are required to fix the 
remaining two coordinates; the condition to determine the foliation of the spacetime by 
a one-parameter family of spacelike hypersurfaces and that to specify the time evolution 
of the radial coordinate. In Eq. ([1]), the radial coordinate has already been fixed so that 
^D~2i^^ '^ becomes the area of the {D — 2)-dimensional sphere labeled by r. Thus, we only 
need the condition for the foliation. As mentioned earlier, we shall consider the foliation 
of this spacetime by the family of maximal hypersurfaces, where a maximal hypersurface 
implies that the trace of the extrinsic curvature vanishes on it, i.e. K = 0. 

In the following, we solve Einstein's equations for the metric ([T]) in the maximal slicing 
condition. From the maximal slicing condition K = 0, we have 

- dt (ln7) + Pl-'dr [In (/327-V2(^-2))] = 0. (2) 

The Hamiltonian and momentum constraints are calculated to give 

{D - = P - 3)(7 - 1) + rdr{ln-f), (3) 

a, [ln(a-i/?7-V^-2^] = 0. (4) 

Using dtK = 0, the evolution equation for the trace part of the extrinsic curvature reduces 
to 

d^a + {D- 2)r~^dra - ^ [5,(ln7)] (9,a = {D - 2)r-^ [{D - 3)(7 - 1) + ra,(ln7)] . (5) 
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The remaining nontrivial component of the evolution equations gives 



dt [In = [{2D - 5)l3-f'^ + {D - 3)a^l3-\^ - 1)] r'^ + 3^-^drf3 

+ \ - 4/57"') 5r(ln7) - (/57"' + a'p-')dr{\na). (6) 

General solutions for a, (3 and 7 are derived as follows. From Eq. (jlj), /3 is determined as 

(3 = T{t)a^r-^''-^\ (7) 

where T(t) is a function of integration. Substituting the above equation into Eq. ([3]), equa- 
tion for 7 is derived to give 



9,(r^- V') = {D- 3)r^"' - p - 1)TV-^, 
and by the integration of this equation, we obtain 

r' = l-h{t)/rf-' + Tyr^^^-^\ 



(8) 



(9) 



where rg{t) is a function of integration. Eliminating (3 in Eq. ([2]) by using Eq. ([7]), and then 
using Eq. ([9]), we have 



(10) 



2T 

Eliminating /3 from Eq. ([6]) and rewriting the result with help of Eq. ffTOl) . we find that is 
a constant: 

9trg = 0. (11) 
Using this fact, Eq. (|TOl) is integrated to give 

dtT r^/" 

Id 

g 

where 



a = /(rg/r;T)V2 



1 + ^ / x^-V(x;T)-^/^rfa; 
'^g ^'^ Jo 



(12) 



fix; T):=l- x""-' + t\^'~''-'^ x'^^'-^l (13) 

If T(t) is determined, 7 and a are subsequently derived by solving Eqs. Qj and f|T2l) . and 
then, /? is determined by Eq. ([7]). Because T(t) is an arbitrary function, we have to impose 
an additional condition for T{t) for specifying a solution. In the next section, we impose 
the refiection symmetry with respect to the wormhole slot and derive the function T{t) that 
specifies this slicing. In Sec. IV, we also consider the case that T(t) is constant and give a 
different class of the maximal hypersurfaces. 



III. REFLECTION SYMMETRIC FOLIATION 



The coordinate system with the choice T = agrees with the Schwarzschild-Tangherhni 
static coordinates, 



1- ^ 



D-3 



dT' + 



l-(^ 

r 



dr"^ + r'^dOi 



2 

D-2- 



(14) 



To derive the coordinates for the general fohation with T 7^ 0, we first prepare the max- 
imally extended Schwarzschild-Tangherlini spacetime (see Appendix A for the Kruskal ex- 
tension), and then, perform coordinate transformation from the Schwarzschild-Tangherlini 
coordinates (r, r) to the coordinates (t, r) for a foliation T 7^ 0. Here, the Schwarzschild- 
Tangherlini time coordinate r is given as r = T{t,r) and satisfies 



dr 
dt 
dr 
dr 



1/2 



.^l/22^^-(I5-2) 



' g ^ 



(15) 
(16) 



Integrating Eq. (fT6|) gives r in the form 

Irg/r 



X 



-dx. 



(17) 



l)/(a;;T)V2 

where X{T) is a function of integration. Substituting the above equation into Eq. ( ITSl) . the 
equation for X is derived as 



dX 

It 



y-i^4-D (x^-3- 1)/(X;T) 



1/2 



,D-3 



dtT 



+ 



X 



f{x-Tfl^ 



dx 



(18) 



Here, we require the hypersurface to have a reflection symmetry with respect to the 
wormhole slot, i.e., r = , which is located in the black hole interior r = rmin < (see the 
Kruskal diagram Fig. [1]); this condition determines T{t). At r = 0, the 1-form V^t, normal 
to the hypersurface labeled by t, should be perpendicular to {d/drY because r is a spacelike 
coordinate in the black hole interior, and thus 
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07 



1/2 







at r = 0. 



(19) 



Because a should be finite everywhere on the slice, the following condition has to be satisfied: 

1 



1 [Train] T) 



0. 



(20) 




maximal slices 



FIG. 1: The reflection symmetric maximal slicing in the Kruskal diagram of Schwarzschild- 
Tangherlini spacetime with D = 5. The maximal hypersurfaces with the reflection symmetry 
with respect to the wormhole slot r = are depicted by solid curves. The dotted curves show the 
r = const, lines, while the dotted straight lines indicate the r = const, lines. The singularity, the 
limit surface, and the horizon are given by r = 0, y/2/2>rg, and rg, respectively. 




Note that the circumferential radius r takes the minimal value r^^^ at r = in the black hole 
interior because of the requirement of the reflection symmetry. Equation (120|) determines 
the value of rmin for a given value of T, i.e., rmin = fminiT). 

The equation /(x; T) = has at most two real positive roots, and the smaller one is 
X = Tg/rmin because f{rg/r;T) = l/7(r;T) holds. By a careful limiting procedure, the 
following fact is found: If X{T) is a smaller root of f{X; T) = 0, then X satisfies Eq. (fTSll ; 
and if X is a larger root, then the integrand in Eq. (fT8l) becomes imaginary that is forbidden 
here. This implies that we have to adopt X(T) = rg/r^i^(T) in Eq. f|T7j) . 

We further require the time coordinate t to agree with the Schwarzschild-Tangherlini time 
coordinate r at spacelike infinity r oo. Then, Eq. f[T7l) gives 



This equation determines the function of T{t). 

If the equation /(x, T) = has a degenerate double root, the integral of Eq. (l2T|l diverges, 
i.e., t = oo. This is the so-called limit surface to which the sequence of maximal hypersurfaces 





(21) 
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asymptotes. Since the root of / = is also the root of the equation df /dx = in this case, 
we obtain 



Urn T\t) = T« 



D-3 
2{D-2) 



D - 1 
_2{D-2) 



(D-l)/(D-3) 



2{D-2) 



The root x = xum of f{x, Too) = is then given by 



2{D 



l/(D-3) 



The minimal radius rmin for t — oo is called the limit radius, and it is given by 



(22) 



(23) 



D - 1 



1 l/iD-3] 



(24) 



2{D-2) 

rn 

In the case of D = 4, this shows the known result r^^ = (3/4)rg [23]. 

Several maximal hypersurfaces with the reflection symmetry at the wormhole slot r = 
for the case of D = 5 are depicted in the Kruskal diagram Fig. [1] (see Appendix A for the 
method of embedding). As in the case D = 4, the sequence of maximal hypersurfaces has 
the singularity avoidance nature for D > 5. 



IV. STATIONARY FOLIATION 



In this section, we turn our attention to a foliation of maximal hypersurfaces which is 
different from the one analyzed in Sec. Ill; foliations for which T(t) is constant and the 
reflection symmetry with respect to the wormhole slot is not imposed in general. As the 
fixed value, we choose T = Too defined in Eq. f l^ . In this case, the sequence of hypersurfaces 
is one-parameter family labeled by t (not by T), and the metric does not depend on the time 
coordinate t. Because the time coordinate basis is not orthogonal to each hypersurface, we 
refer to this foliation as the stationary foliation. 

Several stationary maximal hypersurfaces of T = Too for D = 5 are depicted in the 
Kruskal diagram Fig. [2] (see Appendix A for the method of embedding). The limit surface 
exists also for the stationary foliation with T = Too, and it agrees with the limit surface of 
the refiection-symmetric foliation, studied in the previous section. Figure [2] shows that the 
sequence of maximal hypersurfaces in this class also has the singularity avoidance nature 
and that except for the limit surface, the hypersurfaces of the stationary foliation are not 
refiection symmetric. 
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FIG. 2: The stationary maximal slicings of T = Too in the Kruskal diagram of Schwarzschild- 
Tangherlini spacetime with D = 5 (solid lines). Except for the limit surface r = y^2/3rg, the 
maximal hypersuefaces of T = Too are not symmetric with respect to the wormhole slot. 

The stationary foliation is of special interest in connection with numerical relativity 
performed in the moving puncture approach, because it would be an attractor of the time 



evolution with a dynamical slicing jdjxx = —CaK where C is a constant) 27|] and F-driver 
conditions 2J] as demonstrated in [25] for D = 4 and in 2l| for D = 5. In other words, the 
numerical evolution starting from a hypersurface of stationary foliation has to be unchanged 
in time in these gauge conditions. Therefore, this solution provides a useful benchmark test 
for higher-dimensional numerical relativity codes. 

Black hole simulation in numerical relativity is often performed in the isotropic coor- 
dinates. In Sec. IVA, we describe hypersurfaces of stationary foliation in the isotropic 
coordinates, and study the asymptotic behaviors of the spatial metric. In Sec. IVB, we 
derive the analytic solution of the stationary maximal hypersurface for D = 5 in terms of 
the BSSN variables, and then, demonstrate its usefulness for a benchmark test of numerical 
relativity codes, by performing numerical simulation. 



9 



A. Asymptotic behaviors 



In the isotropic coordinates, the hne element of a (D — l)-dimensional spherically- 
symmetric spacelike hypersurface is written as 

dl^ = ^4/(^-3) + . (25) 

By comparing this with the spatial part of the metric ([T]), the relation between the coordi- 
nates r and R is found: 

^y^dr = V^2/(«-3)^i^^ (26) 

r = ^2/(i>-3)^. (27) 



These equations lead to a differential equation 

dlnR 1 



(28) 



and the formal solution is given by 

i? = i?yM, (29) 

where 

= f . ^ - (30) 

In the following, we analyze asymptotic relations between r and R, and the behavior of the 
conformal factor ip in the distant region r ^ rg and in the neighborhood of r = (given 
by Eq. fl24l) ). one by one. 

In order to study the asymptotic behavior of R in the distant region, we rewrite the 
function /(r) in the form 

/(r) = -ln(rg/r) + F(rg/r), (31) 

where 

F{y) := r dx^—^^^S-. (32) 

Jy x^Jf{x]T^) 

For analyzing its behavior for ?/ — 0, we write F{y) in the form of a Maclaurin series, 

nv) = m - ^(^^ + Oiy'^--'). (33) 

By requiring R to agree with r at spatial infinity r ^ 00, the integration constant Rc is 
chosen to be 

i?^ = rge-^(°). (34) 
10 



Then, the relation between R and r is found to be 

1 



R ^ r 



1 - 



r \ D-3 
I g \ 



for r ^ r„ 



2{D-3) \r 

From Eqs. (!27|) and (!35|) . the asymptotic behavior of the conformal factor ip is given by 

1 /r,\D--i 



(35) 



for i? ^ r„ 



(36) 



which is the well-known relation for D = 4. 

Next, we investigate the asymptotic behavior of R in the neighborhood of r = riim- 
Because /(x; Too) = has the degenerate double root at x = given by Eq. (!23l) . /(x; Too) 
is written in the form 

f{x; Too) = (a^iim - x)'^h{x). (37) 

Here, a function h{x) is a positive definite polynomial of x for < x < xum- Substituting 
Eq. (!37|) into Eq. (!30l) . the function /(r) is rewritten in the following form, 

rfx 1 



Kr) :-- 



where 



■ /r x(xiini x) -i^/7l(xy XUm im J 

i^(l/) : = 



In 



r - rii, 



'^(l - Him/^g) 



Xiii 



Xii 



A//l(Xiim) - Xa//i(x) 



a;(xiim — x)\Jh{x) 



H{rjr), (38) 



(39) 



Note that i?(?/) is finite at ?/ = xum- By taking the second derivative of Eq. ( !37l) . we have 



/i(xiij 



2rfx2 



-(D-l)(D-3) 



2(T'-2) 
T) - 1 



-, (D-5)/(D-3) 



(40) 



Substituting this expression with Eq. (123!) into Eq. (!38|) and using Eq. (l29ll . we obtain the 
relation between i? and r in the neighborhood of r = 



l/y/{D-3){D-2) 



where 



-l/v/(i?-3)(D-2) 



Thus, the asymptotic behavior of the conformal factor is given by 

'nim\(^-3)/2 



R 



for R -C Tg. 



(41) 



(42) 



(43) 
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Although the factor if) becomes steeper near the puncture for higher dimensions, the overall 
conformal factor f[25l) has the universal behavior ~ (rum/i?)^ for 

arbitrary values of D. 

The above result shows that the coordinate origin = has a finite circumferential 
radius r = riim- By contrast, the proper length / from a point labeled with the isotropic 
radial coordinate R to the origin i? = is 

I = r ^ViD-,)^^^^^ ^ ^^^^ p ^ = 00. (44) 
Jo Jo ^ 

This result implies that each stationary maximal hypersurface with T = Too has an infinitely 

long trumpet shape for D > A. 



B. Explicit construction for D = h and numerical evolution 

In the following, we derive solutions of the stationary slice for D = 5 in terms of the 
variables of the BSSN formalism, and then, demonstrate that the solutions are useful for 
checking numerical relativity codes based on the BSSN formalism for D = 5. 

The fundamental variables in the BSSN formalism are different from those in the so-called 



standard 3+1 formalism 28|. In the standard 3+1 formalism, the fundamental quantities 
are the intrinsic metric 'jij, the extrinsic curvature Kij, the lapse function a, and the shift 
vector /?* = 1,...,D — 1 in Z) dimensions). The initial values of these quantities are 
provided by solving the constraint equations, and the subsequent evolution is achieved by 



solving the evolution equations in certain coordinate conditions (see e.g. Ref. [29|]). The 
standard 3+1 formalism prohibits a longterm stable numerical evolution, because constraint 
violation modes grow in the presence of truncation error. In the BSSN formalism, the 
number of dynamical variables is increased to suppress the source of such instability and to 
enable longterm stable simulation. 

The Cgina. fonn of the BSSN fo.,„a..s^ was dese.bed in Refs. fiQ. and it was extended 



to for general dimensionalities in Ref. 



2l[ |. The definitions of dynamical variables in the 



dimensional BSSN formalism are X'- the conformal factor; 74^: the conformal intrinsic 
metric; K: the trace of the extrinsic curvature; Aij: a tracefree extrinsic curvature; and F*: 
a auxiliary D — 1 variable, x jij are defined by 

7ii = X7ii, (45) 
12 



where x is determined so that the determinant of is equal to unity (note that we assume 
to use the Cartesian coordinates). The tracefree extrinsic curvature is defined by 

A^j := X (^K,, - , (46) 

and r* is defined by 

^ (47) 



where 'j^^ is the inverse of i.e., Y ikj = Sj. The variables X; lij-, ^ij-, and are 
evolved, imposing gauge conditions for a and (3\ 

For a numerical relativity simulation, the data of the stationary slice are prepared in the 
Cartesian spatial coordinates {x,y,z,w). In the assumption that the intrinsic metric 'jij is 
conformally fiat, 7jj = 6ij and F* = 0. Furthermore, the maximal slicing condition gives 
K = 0. Since the slice is spherically symmetric, it is sufficient to prepare the data on the 
X-axis, i.e., y = z = w = 0,on which the following relations hold: 

f3^ = py = = (3^ = 0, (48) 

i. . = for z 7^ J. (50) 

Here, R = (a;^ + y"^ + z'^ + w"^)^^"^ indicates the isotropic radial coordinate. This implies that 
the data only for a, Ajui, and x are needed in the spherical polar coordinates. 
For D = 5, it is easy to perform integral in Eq. (1301) to give 

^ - g (3 + ^3 [(rg/r)2 + 3] j (^2{rjr)^ + 15 + 6^2 [{rjr)^ +f]) ' ^^^^ 

The conformal factor x, the lapse function a, and the /^-component of the shift vector are 
given by 

X = ^p-' = {R/rf , (52) 



« = A/l-f^)' + 7^f^)', (53) 



27 V r 
and 

^r. dR 2 rlR 

P = -^T P = — 54) 



The i?i?-component of the tracefree extrinsic curvature is 

2 r 
V3r 



Ann = J. (55) 
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FIG. 3: The values of a, /3^, Ayy, and x on the x-axis for the stationary maximal hypersm'face with 
T = Too (solid curves), and the data after a longterm evolution by the time t = 100 (0), where 
the stationary maximal hypersurface was adopted as the initial data. The units of the length and 
time are rg/2. The data remain approximately stationary after the longterm evolution. 

To describe the data in the isotropic radial coordinate R, we have to give r as a function 
of R. Since the inversion of Eq. (!5T!) cannot be done analytically, we numerically derive the 
relation r = r{R) and then generate the data as functions of R. The analytic initial values 
of a, P^, Ayy, and x on the x-axis are depicted by the solid curves in Fig. [31 Here, we adopt 
rg/2 as the unit of the length. 



We evolve the initial data by using the numerical code recently developed [2l(] and show 
that the data indeed remains stationary in the puncture gauge conditions. We adopt a 
dynamical time slicing condition, 

dta = -2aK. (56) 



This is a simplified version of the 1+log slicing condition 27|] that was studied in Ref. 25 1. 
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FIG. 4: Left panel: The value of ^xx at time t = 60 for resolutions Ax = 0.1 (O), 0.075 (x), and 
0.05 (•). The deviation from the analytic solution 7^"^ = 1 becomes smaller as the resolution is 
increased. Right panel: Snapshots of ^xx for a fixed grid spacing Ax = 0.05 for time t = 20, 40, 
60, and 80. The error gradually increases during the evolution. 



As the spatial gauge coordinates, we adopt the F-driver condition [24 1 



= 2{D -2) ^' and dtB' = dft' - r]B\ (57) 

where t] is constant, chosen to be l/5rg-20/rg. 

The stationary maximal slicing with the isotropic spatial coordinate satisfies the coordi- 
nate conditions ( l56l) and ( 1571) . and thus, the numerical data have to be unchanged during 
numerical evolution with these gauge conditions, if the initial data for x, 7^^, K, Aij, a, 
and /3* agree with those of the stationary maximal hypersurface with T = covered by 
the isotropic coordinate system, together with further initial data 5* = 0. Therefore, the 
analytic solution for these variables given here is used for a test simulation of numerical 
relativity codes based on the BSSN formalism and puncture gauge. 

The data obtained by evolving the above analytic initial data with 5* = is depicted 
by the circles in Fig. [31 In the numerical simulation, the outer boundary is located at 
a sufficiently distant zone (x = 200), and the grid size is uniformly Ax = 0.1. The value 
of 7] is chosen for a wide range as = l/(5rg)- 20/rg, and we confirm that the result does 
not depend on the choice (but we found that numerical error becomes large for very large 
values of if). Figure [3] shows that the numerical data at t = 100 agrees well with the initial 
data t = 0. Thus, we conclude that the numerical data is approximately stationary and 
unchanged in time. 
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The left panel of Fig. H] shows the evolution of numerical values of '^xx on the x-axis 
obtained by r/ = 2/rg for several grid resolutions with Ax = 0.1, 0.075, and 0.05. Because 
the analytic solution is 7^^^ = 1, the deviation from unity indicates the amount of the error. 
This figure clearly shows that the deviation is decreased as the resolution is increased, and 
hence, the deviation from the analytic data is caused only by the numerical error. Here, 
the spatial pattern of the error for a fixed time t with t ^ lOrg was found to depend 
on the grid resolutions. However, as the resolution is increased, the pattern becomes less 
dependent on the resolution and the error decreases approximately at the fourth order. This 
implies that our numerical solution achieves the four-order convergence in the limit Ax 
except at the region near the puncture (where the analyticity of solution is broken and 
the numerical solution should not converge at the fourth order). The right panel shows 
the snapshots of for ^ = 20-80 for the grid size Ax = 0.05. Although the error grows 
during the evolution, the growth rate is small. These results illustrate both the accuracy of 
our numerical simulations and the effectiveness of the foliation by the stationary maximal 
hypersurfaces as a benchmark for checking a five- dimensional numerical relativity code. 



V. SUMMARY 



We have studied the foliation of the D-dimensional Schwarzschild-Tangherlini spacetime 
by the two kinds of one-parameter family of maximal hypersurfaces: the reflection symmetric 
foliation with respect to the wormhole slot and the stationary foliation. We have shown that 
the both foliations have the singularity avoidance nature for D > 5, as in the case of D = 4. 
It is also shown that each hypersurface of the stationary foliation has an infinitely long 
trumpet-like shape in the neighborhood of the black-hole puncture located at the origin of 
the isotropic coordinate. Because the stationary foliation is the attractor of the numerical 
evolution by a dynamical slicing condition [25], both the maximal slicing condition and the 
dynamical slicing conditions will have the preferable nature for the puncture method for 
-D > 4. We presented the explicit solution of the stationary foliation for D = 5, and showed, 
by performing the numerical simulation, that it is useful for a benchmark test of = 5 
numerical relativity codes. 

The remaining issues to be explored are as follows. Although we expect that the puncture 
gauge condition based on a dynamical slicing and F-driver gauge condition also works well 
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for many issues in the higher- dimensional numerical relativity, as demonstrated in Ref. [s^ 
in the four-dimensional case, more detailed studies for a variety of spacetimes are obviously 
needed. For example, it is important to figure out the gauge conditions suitable for simu- 
lating higher-dimensional rotating black hole spacetimes (i.e. Myers- Perry black holes 3l| ) 
and for simulating black holes with high velocity. These are issues to be studied. 
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APPENDIX A: KRUSKAL EXTENSION AND EMBEDDING OF MAXIMAL 
HYPERSURFACES 



In this section, we analyze the Kruskal extension of the D-dimensional Schwarzschild- 
Tangerlini spacetime and explain how to embed maximal hypersurfaces in it. We start from 
the standard metric (fT^ in the static coordinates (r, r) and consider the region r > rg. 
First, we introduce the so-called tortoise coordinate r*. 



dr 



(rg/r) 



D-3 



r -|- 



D-3 



In 



^-1 



G{r) 



Here, defining 



(Al) 



GJr) := cos 



2mT 
D-3 



In 



2r 



cos 



the function G{r) is given by 



2sin 



(D-4)/2 



2mT 
D-3 



2mi 
D-3 

arctan 



cos[2r27r/ {D — 3)] — r/r^ 
sin[2n7r/(D - 3)] 



G{r) = ^2 Gn{r), for even D > 4, 



(A2) 



(A3) 



n=l 
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and 



G{r) = In 



r 

- + 1 



-1 {D-5)/2 

+ J2 ^nir), for odd D>5. 

n=l 



(A4) 



It is easily seen that G{r) is regular for r > 0. Then, we introduce the Kruskal null 
coordinates as 



U = -rgexp[-(D-3)(r-r,)/2rg], 
V = +rg exp [+(/}- 3) (r + r,)/2rg]. 

In these coordinates, the metric ( |T4l) is reduced to the following form: 

2 



(A5) 
(A6) 



ds' 



D-3 



2 -D-3 
^-(D-3)r/rg-G(r) 



J2(^!:^ydudv + r'dn 

n=l 



2 

i:)-2- 



(A7) 



Now, we can extend the spacetime in a similar manner to the four-dimensional case. The 
coordinates U and V introduced by Eqs. (lASP and (]A6P are restricted to the region U < 
and \^ > 0. However, since the metric is regular at ?7 = and V = 0, the spacetime is 
extended to the region U > or V < 0. The maximally extended spacetime consists of four 
regions, and the three regions obtained by the extension are: the black hole region U > 
and V > 0, the white hole region U < and < 0, and the other region U > and < 
outside of the two holes beyond the wormhole slot. The relation between the coordinates 
{U, V) and the static coordinates {t, r) in each extended region is given by appropriately 
changing the sign of Eqs. (lASp and flA6p . Note that the range of r is < r < rg in the 
black and white hole regions while rg < r in the two outside regions. The line element flA7p 
is regular everywhere except at the two physical curvature singularities, r = 0, in the black 
and white holes. See Fig. [2] for the structure of the maximally extended spacetime. 

In order to embed maximal hypersurfaces in the Kruskal diagram, we consider the coor- 
dinate transformation between the Kruskal coordinates {U, V) and the coordinates {t, r) for 
maximal slicing. Since {U, V) and (r, r) are related through Eqs. flASp and flA6p whereas 
(r, r) and (t,r) are related through Eqs. ffT^ and (ITBl) . we obtain 

-1 



dU 
dr 

dV 

dr 



D-3 


1 - 




2r 

t ^'g 






D-3 


1 - 




2r 







(i + y/2r^-(^-2))f/^ 



(A8) 
(A9) 
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Once the arbitrary function of integration T{t) is determined, we obtain a maximal hyper- 
surface in the Kruskal diagram by integrating these two equations. 



[1] F. Pretorius, Phys. Rev. Lett. 95, 121101 (2005) |arXiv:gr-qc/05 07014]. 

[2] M. Campanelli, C. O. Lousto, P. Marronetti and Y. Zlochower, Phys. Rev. Lett. 96, 111101 

(2006) [ arXiv:gr-qc/0511048| . 
[3] J. G. Baker, J. Centrella, D. I. Choi, M. Koppitz and J. van Meter, Phys. Rev. Lett. 96, 



111102 (2006) |a rXiv:gr-qc/0511103] . 
[4] F. Herrmann, D. Shoemaker and P. Laguna, |arXiv:gr -qc/0601026[ 
[5] U. Sperhake, Phys. Rev. D 76, 104015 (2007) [arXiv:gr-qc/0606079| . 

[6] M. A. Scheeh H. P. PfeifFer, L. Lindblom, L. E. Kidder, O. Rinne and S. A. Teukolsky Phys. 

Rev. D 74, 104006 (2006) |a rXiv:gr-qc/06070 56|. 
[7] B. Briigmann, J. A. Gonzalez, M. Hannam, S. Husa, U. Sperhake and W. Tichy, Phys. Rev. 

D 77, 024027 (2008) |arXiv:gr-qc/0610128l . 
[8] M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 (1995). 

[9] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 024007 (1998) [arXiv:gr-qc/9810065| . 
[10] T. Banks and W. Fischler, |arXiv:hep-th/9906038 l 

[11] S. Dimopoulos and G. Landsberg, Phys. Rev. Lett. 87, 161602 (2001) [arXiv:hep-ph/0106295] . 



[12] S. B. Giddings and S. Thomas, Phys. Rev. D 65, 056010 (2002) [arXiv:hep-ph/0106219| . 
[13] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 
|arXiv:hep-ph/98033"l5] ; 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvah, ibid. 436, 257 (1998) 
|arXiv:hep-ph/98043"98] . 



[14] L. Randah and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) |arXiv:hep-ph /9905221]; 
[15] P. Kanti, Lect. Notes Phys. 769, 387 (2009) jarXiv:0802.2218l [hep-th]]. 

[16] U. Sperhake, V. Cardoso, F. Pretorius, E. Berti and J. A. Gonzalez, Phys. Rev. Lett. 101, 

161101 (2008) [a rXiv:0806. 17381 [gr-Qc]]. 
[17] M. Shibata, H. Okawa and T. Yamamoto, Phys. Rev. D 78, 101501(R) (2008) [ arXiv:0810.4735] 

[gr-qc]]. 

[18] U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, T. Hinderer and N. Yunes, arXiv: 09071252] 



19 



[gr-qc]. 

[19] H. Yoshino and Y. Nambu, Phys. Rev. D 67, 024009 (2003) |ar Xiv:gr-qc/0209 003]. 
[20] H. Yoshino and V. S. Rychkov, Phys. Rev. D 71, 104028 (2005) |arXiv:hep-th/050317T] . 
[21] H. Yoshino and M. Shibata, arXiv:0907.2760 [gr-qc]. 
[22] F. R. Tangherhni, Nuovo Cim. 27, 636 (1963). 

[23] F. Estabrook, H. Wahlquist, S. Christensen, B. DeWitt, L. Smarr and E. Tsiang, Phys. Rev. 
D 7, 2814 (1973). 

[24] M. Alcubierre, B. Briigmann, P. Diener, M. Koppitz, D. Pollney, E. Seidel and R. Takahashi, 

Phys. Rev. D 67, 084023 (2003) |arXiv:gr-qc/0206072] . 
[25] M. Hannam, S. Husa, B. Briigmann, J. A. Gonzalez, U. Sperhake and N. O. Murchadha, J. 

Phys. Conf. Ser. 66, 012047 (2007) [ arXiv:gr-qc/061 2097|. 
[26] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 024007 (1998) [arXiv:gr-qc /9810065|. 
[27] C. Bona, J. Masso, E. Seidel and J. Stela, Phys. Rev. D 56, 3405 (1997) |arXiv:gr-qc/97090T6] . 
[28] J. W. York, Jr., in Sources of Gravitational Radiation editted by L. L. Smarr (Cambridge 

University Press, 1979), pp. 83. 
[29] L. Smarr and J. W. York, Jr., Phys. Rev. D 17, 2529 (1978). 

[30] M. Hannam, S. Husa, D. Pollney, B. Briigmann and N. 6 Mmxhadha, Phys. Rev. Lett. 99, 

241102 (2007) [a rXiv:gr-qc/0606099 |. 
[31] R. C. Myers and M. J. Perry, Annals Phys. 172, 304 (1986). 



20 



